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We calculate the low- frequency optical conductivity o~(lo) of clean metals and superconductors 
at zero temperature neglecting the effects of impurities and phonons. In general, the frequency 
and temperature dependences of a have very little in common. For small Fermi surfaces in three 
dimensions (but not in 2D) we find for example that Recr(o; > 0) » const, which corresponds to a 
scattering rate F oc lo 2 even in the absence of Umklapp scattering when there is no T 2 contribution 
to r. In the main part of the paper we discuss in detail the optical conductivity of d-wave super- 
conductors in 2D where Kea(io > 0) oc to 4 for the smallest frequencies and the Umklapp processes 
typically set in smoothly above a finite threshold loq smaller than twice the maximal gap A. In cases 
where the nodes are located at (±7r/2, ±7r/2), such that direct Umklapp scattering among them is 
possible, one obtains Rea(ui) oc lo 2 . 

PACS numbers: 



Introduction: Optical conductivity is a powerful tooli 
to study the properties of a strongly correlated metal. 
The frequency dependence in particular can give detailed 
information on the excitation spectrum of a system (gaps, 
phonons, magnons, interband transitions, ...) which in 
general cannot be extracted from, for example, the tem- 
perature dependence of the conductivity. 

In a superconductor the electronic contribution to the 
optical conductivity Re o~(lo) can be separated - at least 
in simple situations - into three different contributions. 
First - and most important - superconductivity implies 
the existence of a (5-peak at lo = whose (Drude-) weight 
is given by the condensate fraction. Second, thermal ex- 
citations at small but finite temperatures, T > 0, are ex- 
pected to lead to a rather sharp peak centered again at 
u> = 0, whose width is identified with the scattering rate 
of the thermal excitations and is strongly temperature 
dependent (in cases where impurity scattering can be 
neglected). Finally, all other contributions at finite fre- 
quency are usually called "incoherent background" . This 
background depends only weakly on temperature T. It is 
the goal of this paper to discuss the low-frequency prop- 
erties of this incoherent background. More precisely, we 
consider the optical conductivity for frequencies lo > at 
T = when thermal excitations are absent. 

In most experimentally relevant situations the opti- 
cal conductivity of (conventional) superconductors at low 
frequencies is dominated by elastic impurity scattering. 
The theory of optical conductivity in such systems was 
developed very early by Mattis and Bardcen 2 . Inelastic 
scattering is more important in strongly interacting su- 
perconductors and rather clean samples, and therefore 
the optical conductivity in d-wave superconductors has 
been studied quite extensively in th e con text of high-T c 
superconductors (see for example [il l llo| and references 
therein). Motivated by experiments, these investigations 



have mainly investigated the influence of scattering from 
collective modes and the interplay with impurity scat- 
tering. In this paper we systematically investigate the 
zero-temperature optical conductivity arising from the 
interplay of band-structure effects and electron-electron 
interactions taking into account all relevant vertex cor- 
rections. While the main focus of this paper is the investi- 
gation of d-wave superconductors, we also briefly discuss 
the optical conductivity of clean Fermi liquids and s-wave 
superconductors. 

Method: According to the Kubo formula, the optical 
conductivity is given by 
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where ((J, J)) w is the the current-current correlator, 
«J, J)) w = -i/ oo dte i ( w + i °) t ([J(t), J(0)]). When calcu- 
lating the optical conductivity perturbatively, it is im- 
portant to take into account both vertex and self-energy 
corrections. For example, in a Galilean invariant system 
with a quadratic dispersion, e k = fc 2 /(2rn), vertex and 
self-energy corrections cancel exactly as the total electri- 
cal current is a conserved quantity. But even in clean 
non-Galilean invariant systems, i.e. for electrons mov- 
ing in a periodic crystalline potential, massive cancella- 
tions between self-energy and vertex corrections occur, 
especially if there is little Umklapp scattering close to 
the Fermi surface. To take into account vertex and self- 
energy corrections on the same footing, one in general 
has to solve an integral equation (a vertex equation or, 
equivalently, a linearized quantum Boltzmann equation) 
to obtain the correct conductivity even to lowest order in 
perturbation theory. 

However, at zero temperature and in the absence of 
disorder one can avoid the substantial technical difficul- 
ties involved in solving multi-dimensional integral equa- 
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tions by the following argument: In general, one can 
express the conductivity in the form <j{lj) = T (J)-iu> 
where \ 1S identified with the total optical weight and 
the (frequency-dependent) scattering rate Re T(u>) can be 
calculated from the integral equations described above. 
However, for <C lu this simplifies after multiplica- 

tion with lu 2 to 

cu 2 Re a(u) = lu 2 Re — « x Re T(w) (2) 

I (w) — ilu 

Note that there is no contribution from the 5-function at 
uj = due to the ui 2 prefactor. For weak interactions T is 
small and therefore we can obtain o~(lo > 0) in a straight- 
forward perturbative expansion, i.e. without solving any 
integral equations, from the right-hand side of 

Im id.J,d f J\ 

Re a(u > 0) = VW I { 11 " , (3) 

or 

provided that |r(w)| <C cu. As d t J is already linear in 
the interactions (see below), it is sufficient to leading or- 
der to evaluate the correlation function in © to zeroth 
order in the couplings. We will use this approximation 
only at T = 0. At any finite temperature, the scatter- 
ing rate Y{uj — * 0) is constant and therefore the method 
described above will break down for lu —* but remains 
valid at higher frequencies where |r(w)| <C u>. Note that 
within the so-called memory-function approach^ one uses 
essentially identical formulas to calculate T(lu). 

If one is only interested in the qualitative behavior of 
Recr(a>) at low u, i.e. in the power-law obtained for 
lij — > 0, one can relax the condition |r(u;)| <C lu and 
replace it by Rer(w) <C (1 — c)u> for lu < lu* where lu* 
is some characteristic frequency and c (which can be of 
order 1) is obtained from IuiT(lu) « clu for lu <ti w* . 
As the latter condition is fulfilled in all cases discussed 
below, we expect that all our results are qualitatively cor- 
rect at sufficiently low frequencies (a possible exception is 
discussed below) even in a strongly interacting systems. 
(Backflow and other Fermi liquid renormalization effects 7 
will only change prefactors, and multi-particle scattering 
processes are suppressed for w->0 due to the restricted 
phase space.) 

In the following, we will first consider the optical con- 
ductivity of a clean Fermi liquid at T — 0. This will serve 
as a reference for our results on d-wave superconductors 
presented in the second part. 

Metals: In a one-band model, the electrical current is 

given by J = J2ka v k c Lx c kcr ^ n * nc P resence 01 interac- 
tions and in the absence of Galilei invariance the current 
is not conserved with d t J = iEkk'q^' ^ q ( v k + v k' - 

v k+q- v k'- q ) c L c k+ q , < x c L ff ' c k'-q, ( T' for a density-density 
interaction U . The change of current is proportional to 
the difference of incoming and outgoing velocities. In the 
following we will assume that the (screened) interaction 
U ~ U depends only weakly on the transferred momen- 
tum q. 



For a (normal) metal we therefore obtain at low fre- 
quencies and T = using Eq. © 

Re a{uJ > o) « ±J! ]T / x / 3 (l - / 3 )(1 - / 4 ) 

1234G 

X « + V 3 ~ V 2 ~ O 2 < 5 l+2,3+4+G 

x [5(u-{e A + e 3 -E 2 -e 1 )) - (w -u)] (4) 

where 1,...,4 denote the momenta k 1 ,...,k 4 in the first 
Brillouin zone and momentum is conserved modulo re- 
ciprocal lattice vectors G. To perform the momentum 
integrals it is useful to split into a component perpen- 
dicular to the Fermi surface and an angular integration 
parallel to it. For small to only a thin shell of width to jv F 
contributes for each of the three relevant momentum in- 
tegrations, implying an to 3 dependence which cancels the 
1/ui 3 prefactor. 

If the Fermi surface is sufficiently large such that Umk- 
lapp scattering processes can take place, one therefore 
obtains the well-known result that 

Re cr(ui > 0) « const., r(w) oc ui 2 (5) 

as the 4 velocities sum up to a finite value of the or- 
der of v F in this case. The constant "incoherent back- 
ground" corresponds according to Eq. @ to a scatter- 
ing rate T(u>) oc lu 2 characteristic of a Fermi liquid with 
Umklapp scattering in 2 or 3 dimensions. 

Less well known is the corresponding result for a small 
Fermi surface (k F < G/4) where Umklapp scattering at 
the Fermi surface is not possible. Here the situations in 
2 and 3 dimensions are quite different. For a generic (not 
too complex) Fermi surface in 2D, momentum conserva- 
tion in the limit lu — > can only be fulfilled by choosing 
kj^ = — k 2 and k 3 = — k 4 (or kj^ = k 3 ^ 4 and k 2 = k 4 ^ 3 ). 
Therefore the sum of the velocities also vanishes linearly 
in lu for lu — > and one obtains from power counting 

Rc <t(lu > 0) cx lu 2 , T(oj)(xlu 4 (6) 

for a small Fermi surface in d — 2. 

The situation is quite different for a system with a 
small Fermi surface in 3D, where momentum conserva- 
tion on the Fermi surface does not require that the rele- 
vant moments are located opposite to each other. There- 
fore the sum of the 4 velocities in Eq. wm generically 
not vanish and one finds 

Re o(lu > 0) w const., r(w) cx lu 2 (7) 

for a small Fermi surface in d = 3: Even without Umk- 
lapp processes the scattering rate varies as T(u) oc lu 2 \ 
Note that the frequency and temperature dependence 
of the conductivity are drastically different in this case. 
T(T) does not vary as T 2 but the two-particle scattering 
rate is exponentially suppressed; multi-particle processes 
lead to a power law Y oc T 2 ™~ 2 where the integer n de- 
pends on the size of the Fermi surface^, n ~ G/(2k F ). 
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The disparate behavior of T(u,T = 0) ~ U 2 uj 2 and 
r(w = 0, T) ~ J7 2 e~ A / T + u n T 2n - 2 can easily be under- 
stood once one realizes that, in the absence of Umklapp 
scattering, on the one hand the current is not conserved 
while on the other hand the momentum is conserved. 
As explained in detail e.g. in Refs. 00 > the component 
of the current "perpendicular" to the momentum does 
decay rapidly giving rise to the frequency independent 
incoherent background of Eq. J7J. The dc conductiv- 
ity is, however, determined by the long-time decay of 
the component of the current "parallel" to the momen- 
tum and therefore by the decay rate of the momentum, 
i.e. by Umklapp processes which are very rare for small 
Fermi surfaces. It is likely that the rather general results 
Eq. H6I7(I have been discussed before in the literature but 
we are not aware of a directly relevant reference. 

It should be clear from the discussion of Eq. (|7J) given 
above that also in the presence of Umklapp scattering, 
when r(w,T) ss a(k B T) 2 + b(hw) 2 there is in general no 
simple relation between the constants a and b. We em- 
phasize this fact as in the experimental literature such 
a relation has sometimes been claimed to exisfei*iS*ii but 
is actually not observedii*i£. Note that recent progress 
in the experimental methods allows precise measure- 
ments of the optical conductivity at low frequencies and 
temperatures^. 

Superconductors: We now turn to the calculation of the 



T = optical conductivity in superconductors neglecting 
again phonons and impurities. Our main interest is the 
case of a d-wave superconductor in d = 2 on a square 
lattice with unit lattice spacing. To describe the super- 
conducting state we use weakly-interacting Bogoliubov 
quasiparticles (QPs), d kcr = u k c\ a - crv k c_ k -, which di- 
agonalize the BCS Hamiltonian H BCS = J^ka £ k c ka c ka~\~ 

E k A k4 T c -fci + h - c - = E kCT E k d L d k a wherc C L is thc 
electron creation operator and E k = \J e\ + A| the BCS 
energy. The electric current is given by 

J = E W k C la C ka = E V k d L d ka » ( 8 ) 
kcr ka 

where it is important to realize that it is the bare veloc- 
ities v k = de k /dk rather than dE^/dk that enter if the 
current is expressed in terms of the BCS quasiparticles. 

Within the BCS approximation the current is con- 
served, [J, H BCS ] — 0, and there is no optical weight at 
finite frequencies. To calculate the optical conductivity 
it is therefore essential to include the interaction of the 
quasiparticles. The Hamiltonian for the QPs is given by 
H = H BCS + H lnt , where the (properly normal ordered) 
local density-density interaction H int = 2U'^2 i n i ^n i ^ can 
be rewritten as 



^int = U J2 r l3 r 24 d h d il d ll d ^ + h - C - + 2f 12 r 34 d h d ll d L d la + h ' C - + r l2 T 34 d l<J d \ S d 2a d lo + f 14 f 23 d L d L' d 2a' d la 

(9) 

where = = %i> 2 + w i u 2' ^ij ~ ^ji ~ u i u 2 ~ v i v 2 an< ^ * = The momentum sums conserve crystal momentum 
and the spin sums are only over repeated indices. This expression can be derived by keeping the fluctuations around 
mean-field theory in the BCS approach. The various terms describe not only the scattering of quasi-particles (and 
holes) but also the breaking up and recombination of Cooper pairs. 

While [J, if B cs] = 0> * ne curren t J decays in the presence of the interactions between QPs, d t J — — i[J, H int ]. It is 
now straightforward (albeit somewhat tedious) to evaluate the contributions to to lowest order in the interactions, 
and we obtain 

Re <r( u ) = ^ + <^ q M + <>) - (u, <- -«)) (10) 

= E ( r 12 r 34 " r 13 r 2if S l+2+3+4-G K + V 2 + V 3 + V %f 
1234G 

x [(1 - /,)(1 - / 2 )(1 - / 3 )(1 - /J - fJ 2 f 3 f A ] 5(uj - (E 1 +E 2 + E 3 + EJ) 

Cq( W ) = 4 E ( f 12 r 34 ~ f 14 r 2 3 ) 2 K ~ V 2 ~ V 3 ~ V lf S l-2-3-4+G 
1234G 

x [/,(1 - / 2 )(1 - / 3 )(1 - /J - (1 - /J/a/g/J 6(uj - (-E, +E 2 + E 3 + EJ) 

<qH = E (1 " - h)hh K + «? - 4 - O 2 ^l + 2-3-4 +G 

1234G 

X [ 4 ( r i 2 r 34 + f 14 f 2 3 ) 2 + 2(^14^3 ~ Wllf] 5 (<" ~ (^1 + E 2 ~ E 3 ~ E i)) 

I 



The first (second, third) contribution comes from the first (second, last two) scattering terms in the Hamil- 
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tonian JJjJ. At zero temperature, obviously only the first 
term <fi (oj) survives as all QPs have positive energies 
and/ i = /( J B ki ) = 0atT = 0. 

In the case of an ultra-clean s-wave superconductor, 
a direct consequence of Eq. (|10fl is that the gap in the 
optical conductivity is of size 4A while it is 2A for dirty 
superconductors^ as has previously been noted by Oren- 
stein et al£ Obviously one has to ask whether this re- 
sult will also hold to higher order in perturbation theory. 
To answer this question, one has to investigate whether 
symmetries and corresponding selection rules allow for 
an optical transition from the ground state of the su- 
perconductor to a 2-quasiparticle excited state by an 
operator of the form J^kk' a kk'^kcr^-k'(r' where t? are 
the creation operators of the fully renormalized "true" 
quasiparticles of the system (which can only be identi- 
fied with the BCS quasiparticles for weak interactions). 
Symmetries strongly restrict the form of ■ Transla- 
tional invariance on the lattice, for example, implies that 
a kk' = a k CT ^(k — k') in the absence of impurities with 
a™ — — as the quasiparticles are fermions. If the 
superconductor does not break time-reversal invariance 

one has (a^)* — ) ancl ( a k )* = ~( a u ) ano ^ m a 
crystal with inversion symmetry one has a^ 7 = a£ " . 
In the absence of spin-orbit coupling, i.e. if spins are 
rotationally invariant one finds that ajj = — and 
a™ = 0. From this we can conclude that, in the absence 
of disorder and in the presence of inversion symmetry, 
a™ vanishes and the optical gap is therefore 4A for an 
s-wave superconductor in the absence of spin-orbit cou- 
pling. In the presence of impurities, however, the gap2> 
is only 2A. Interestingly, the symmetry analysis sug- 
gests that even in a generic inversion-symmetric clean 
crystal, high-order spin-orbit processes could possibly in- 
duce relevant low-energy processes not included in 11U|) 
which lead to a gap of size 2A. All the low-order re- 
sults presented below may therefore not be valid in the 
presence of sizable spin-orbit coupling. Note also that 
phonons and other low-energy collective modes with en- 
ergies smaller than 2A can induce optical weight in the 
frequency window 2A < u> < 4A. The precise functional 
form of the optical conductivity of an s-wave supercon- 
ductor for oj > 4A will not be discussed in detail here. 
It depends on the dimension and on the angular depen- 
dence of A. Generically the onset will be smooth and 
of the form (oj — 4A) 2 . Therefore a precise experimen- 
tal determination of A using a feature close to 4A will 
be rather difficult. For all conventional s-wave supercon- 
ductors we anyhow expect that impurity scattering will 
dominate even for the cleanest available samples leading 
to the well-known 2A gap which is much easier to detect. 

A d-wave superconductor in two dimensions (as re- 
alized in high-temperature superconductors) with point 
nodes along the diagonals of the quadratic Brillouin zone 
has a vanishing gap in nodal direction. For frequencies 
small compared to the maximal gap A, oj <C A, all QPs 
are created in the vicinity of the nodes, so we expand 




FIG. 1: The four possible types of scattering process at T = 
0, in which four QPs are created. The circle represents the 
Fermi surface, the solid arrows the QP momenta, and the 
dashed arrows the reciprocal lattice vector G. The shaded 
region indicates the size of the superconducting gap and the 
ellipses a constant energy contour in the vicinity of each node. 

the dispersion around them. Writing k = k d + k the 
most generic band structure consistent with the square 
symmetry of the lattice is 

v 1 

+ + «*) + Fk x k v {k x + Ky ) + 0( K 4 ) (11) 

where the constants D, L, F determine the deviation of 
the dispersion from that of a free-electron gas and m* is 
an effective mass. 

There are four qualitatively different terms that ap- 
pear in the sum for <p'^ p in l)l()[l. which are sketched in 
Fig. (i) G = (27r, 2ir) and hence all four QPs in one 
node; (ii) G = (2tt, 0) and two QPs in each of two "per- 
pendicular" nodes; (hi) G = and one QP in each node; 
(iv) G — and two QPs in each of two opposite nodes. 
These give rise to very different dependences on oj and 
doping, as we now discuss. 

The role of Umklapp scattering is determined by the 
distance of the nodes from (ir/2,ir/2) which we denote 
by 
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In processes of type (i) the four QPs have very simi- 
lar velocities (recall that it is the normal-state velocity 
that contributes) and so the contribution to <p''(u>) is 
large. However, it is only possible to create four QPs of 
arbitrarily low energy if the nodes are situated exactly 
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Fermi surfaces (e.g. optimally doped Bi-2212 according 
to Ref. 0, see also Fig. |3J), the gap for Umklapp pro- 
cesses in a d-wave superconductor will be smaller than 
2A, 



lu q < 2A 



(15) 



as typically an Umklapp process will exist where two QPs 
are located at the nodes and the two other somewhere 
else on the Fermi surface. 

To obtain the frequency dependence, we ignore in a 
first step the coherence factors and velocity prefactors 
and evaluate the integral 



FIG. 2: Zero-temperature optical conductivity, Kea(cj > 0), 
for an ultra-clean d-wave superconductor from a numeri- 
cal evaluation of Eq. (1 10> for a model with nodes close to 
(7r/2,7r/2) [A = 0.31, v F = l,v F /v A = 5, ui = 0.1 and 
lu'o = 0.02]. Inset: logarithmic plot of Rea(uj > 0)/oA The 
thin lines show the contributions from the various processes 
shown in Fig. Processes of type (i) and (ii) are gapped by 
u>o and cu' , respectively. Due to the smooth onset of type-(i) 
Umklapp processes, Eq. 12011 . there is almost no feature at t^o- 
The numerical results are fully consistent with the power-laws 
of Eqs. G7HHJ. 



at (7r/2,7r/2); otherwise there is an excess momentum 
4<5fc nodo which must be carried by the QPs, so that at 
least one of them is situated a finite distance from the 
node. Accordingly, absorption can only occur for fre- 
quencies above the threshold 



Av F Sk n 



)de ' 



(13) 



Processes of type (ii) resemble those of type (i) , since 
again the velocities add. However, the fact that the nodes 
are at right angles to one another reduces the threshold 
frequency as the excess momentum ( fc noc i ) ~ G = 
(2V2Sk nodc ,0) = V26k node [(l,l) + (T-l)] can be split 
into two components parallel to the Fermi surface at the 
nodes where the velocity of the QPs v A — dE k /dk^ = 
dA^/dfcy is much smaller. This leads to a considerably 
smaller threshold frequency 



u Q m <±v A 6k node 



(14) 



J2 5 (w - {E l + E 2 + E 3 + E_ ( 1+2+3) ) ) 

1,2,3 



(^F (16) 



for Sk nodc = and small to in each of the 4 cases shown 
in Fig. ^ This can be done by scaling the momenta 
perpendicular and parallel to the Fermi surface at the 
node by lo/v f and u>/v A , respectively. In the cases (i) 
and (iv) shown in Fig. when all nodes are parallel 
to each other, c, and c 4 are constants of order 1. The 
situation is slightly more complicated in the cases (ii) and 
(iii) where by choosing a proper rescaling procedure we 
find c 2 ~ c 3 ~ v A /v F . 

Eq. (|16|l does not include the effect of the velocity pref- 
actor (J2 v i) 2 anc > of the combination (r, 
of coherence factors in Eq. (|10fl . At the nodes, the 

herence factors y'l ± are rapidly varying functions 

of order 1. For 8k„ 



13' 24' 



co- 



-nodo = t ne y change the result only 
quantitatively but not qualitatively (as we have checked 
numerically) but can become important for Sk nodc =/= 
as discussed below. For the Umklapp processes (i) and 
(ii) the velocities just add up to give a finite prefactor 
of order vp. If the nodes are located at (tt/2, 7t/2), we 
therefore obtain 



<j(lo) (X 



U 2 



J F A 



type (i) for 5k 



node 



(17) 



and similarly 



cr(uj) cx 



U 2 



F A 



type (ii) for 5k n 



= 0. 



(18) 



where this simplified formula is only valid if Sk nodc < These power- laws can also be observed for to > uj ,oj' if 
A/w F when corrections to the Dirac spectrum close to the 
nodes can be neglected. Note that the construction de- 
scribed above reduces the available phase space for scat- 
tering and so the contribution close to uj' is smaller than 
that of type (i) processes by a factor of v A /v F . 

In most realistic situations (including most of the 
cuprates) the point node will not be located close to 
(7r/2,7r/2) and ^fc nodc will be larger than A/u„. In 
this case the gap for Umklapp processes will depend 
on details of the band-structure. For sufficiently large 



Afc nodo is finite but small. This can be seen in the inset 
of Fig. [2J which discusses the various regimes based on a 
numerical evaluation of Eq. (|1C)|> . 

Due to momentum conservation, the leading contribu- 
tion to ^2 vanishes for the non-Umklapp processes (iii) 
and (iv). But band-structure effects break Galilean in- 
variance and one obtains a low-energy contribution even 
in the absence of Umklapp. The leading term is given 

Dn y . Although 



by Dkk, in (tTTt which leads to w 

x y 



the sum J^. still vanishes at this order for processes of 
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coTeVl 

FIG. 3: Optical conductivity, Kea(ui > 0), for an ultra-clean 
d-wave superconductor (dashed line) using the bandstructure 
of optimally doped Bi-2212 taken from |l3j and a d-wave gap 
of size A = 0.028 eV (arbitrary units on the y-axis). For ref- 
erence, the solid line shows the T = optical conductivity in 
the normal state (A = 0) which is constant for low frequencies 
due to Umklapp processes, see Eq. |0. The 5-peak at uj — is 
not shown. For lowest frequencies, uj < loo, one finds a oc to 4 
in the superconducting phase as can be seen on the logarith- 
mic scale of the inset. This regime is, however, practically 
not observable due to the small prefactor. Instead one finds a 
very smooth onset [see Eq. 1211 1 for 10 > u>o * 0.05 eV < 2A 
I15H . Note that the nodes are not close to (7r/2,7r/2) for the 
bandstructure considered in this figure. 

type (iv) , it remains finite if the geometry is determined 
by (iii) and we obtain from a scaling analysis 

U 2 D 2 

a(uj) oc — — g- w 4 type (iii) for uj — > 0. (19) 

V F V A 

2 

While this term is suppressed by the tiny factor ^ ^ 

compared to l|17fl . it is nevertheless the leading uj — > 
correction when the nodes are located away from 
(7r/2,7r/2). Eq. (|19|l therefore describes the typical low- 
frequency optical conductivity of a 2-dimensional d-wave 
superconductor in the absence of impurities (c.f. insets 
of Figs. 1213)) . Processes from the scattering geometry (iv) 
are always subleading and only give rise to contributions 
oc uj 6 . It is worth noting that the prefactor of ltl^|) - 
not shown in the equation - turns out to be numerically 
very small, approximately a factor of 20 smaller than the 
prefactor of l|17|) and more than a factor of 100 smaller 
than the corresponding numerical prefactor of 1(18(1 if we 
assume a local interaction U. In general completely dif- 
ferent matrix elements enter the various scattering pro- 
cesses (i) - (iv) and therefore their relative magnitude 
depends on details of the relevant interactions. But the 
smallness of the contribution may imply that in actual 
measurements the low-frequency uj 4 regime is never ob- 
servable, see Figs. |21 and 

As the non-Umklapp contribution l|19fl to the opti- 



cal conductivity is very small and difficult to detect 
experimentally, it is worthwhile to investigate the pre- 
cise form of the onset of Umklapp terms at uj > uj q , uj' q 
in the generic case when the nodes are not located at 
(7r/2,7r/2). Consider for example the scattering geom- 
etry (i) in Fig. ^ At u> — u> the components Kn 
of all 4 momenta parallel to the Fermi surface will be 
zero, so e k = E, and therefore the coherence factors 
( r i2 r 34 ~ r i3 r 24) °f ^1- P^ft wu l vanish. As a con- 
sequence the onset of Umklapp processes will be very 
smooth and of the form 

U 2 

cr(uj) oc 5- (uj — uj n ) 2 , type (i) for uj > uj n (20) 

v F v% 

and 

U 2 

<t(uj) oc 2 2 (oj — lj' q ) 2 , type (ii) for oj > uj' q (21) 
v F v A 

as we have checked numerically, see Figs.|3]and|2] Formu- 
las for uj q and oj' are given in Eqs. ii 131 115(1 above. The 
prefactors in Eqs. 1(20121(1 are only valid for very small 
uj 01 uj' <C fcpU^/fp when one can use a Dirac spectrum 
for the nodal quasiparticles; however, the frequency de- 
pendence close to the onset frequency is also quadratic 
for larger values of uj and oj' as we have again checked 
numerically for example in Fig. [3] which shows the optical 
conductivity in a model which uses the bandstructure 13 
of Bi-2212. 

All results shown above rely on the fact that at lowest 
energies the nodal dispersion takes the form of a Dirac 

cone, E k = ^J(v F k ± ) 2 + (v A k^) 2 . But already at a very 

low energy scale, E c = m*v A /2 ~ A 2 /e F <§; A, one 
has to take into account the curvature of the Fermi sur- 
face which bends contours of equal energy into a ba- 
nana shape. It is therefore important to check which 
of the results calculated above remain unmodified at this 
crossover scale - the existence of such a small energy 
scale will otherwise make the experimental determina- 
tion of power laws extremely difficult. Fortunately, it 
turns out that our results in the scattering geometry (ii) 
and (iii), i.e. Eqs. 1(18119(1 . are not affected by E and re- 
main valid up to energies of the order of the maximal gap 
A. This can most easily be seen by rewriting momentum 
conservation in polar coordinates while scaling k — k F 
with Lu/v F and the polar angle </> with oj/(k F v A ). Using 
the same analysis for geometry (i), one finds a crossover 
at the energy E c and Eq. I(17|) has to be multiplied by a 
factor EJuj for E c < w < A. 

Conclusions: The frequency dependence of the opti- 
cal conductivity at zero temperature and finite frequen- 
cies describes how the electrical current can decay. The 
example of a Fermi liquid with a small Fermi surface 
shows that the temperature and frequency dependencies 
of (j(u>, T) have very little in common and may result from 
completely different processes. The zero-temperature op- 
tical conductivity of d-wave superconductors turns out 
to be rather complex even for frequencies much smaller 
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than the maximal gap. While we hope that our calcu- 
lation can serve as a reference for the interpretation of 
the incoherent background, a direct observation of the 
predicted power-laws will be difficult as the calculated 
contributions turn out to be both small in size and very 
smooth in their frequency dependence (see Fig.P|and|^)l. 
Therefore it will be very difficult even in very clean crys- 
tals to separate the predicted effects from the effects of 
elastic impurity scattering. 

An interesting open question is whether and how spin- 
orbit interactions modify the results presented in this pa- 
per. Based on a symmetry analysis, we argued that spin- 
orbit interactions can open new channels for current re- 
laxation in a superconductor - even in the presence of in- 
version symmetry. Neglecting such relativistic effects, we 
believe that our results are valid even in strongly interact- 
ing superconductors at sufficiently low frequencies when 
multi-particle scattering is suppressed due to phase-space 



restrictions. This will also be the case if the interactions 
are mediated e.g. by (short-ranged) spin-fluctuations^, 
provided the system is not located directly at a quantum- 
critical point. 

At small but finite temperatures thermal excitations 
induce a characteristic sharp peak in the low-frequency 
optical conductivity. The calculation of this prominent 
feature in a d-wave superconductor taking into account 
the relevant vertex corrections and approximate conser- 
vation laws&*i is left as a challenge for the future - while 
the T — results of this paper can provide a reference for 
this calculation, the simple methods used here will not 
be sufficient to describe the finite temperature regime. 
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